Unusual ferromagnetic superexchange in CdVC^: The role of Cd 
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A microscopic magnetic model of the low-pressure modification of CdVOs is established, based 
on density functional theory (DFT) band-structure calculations, magnetization measurements, and 
quantum Monte-Carlo simulations. This compound is a rare example of a quasi-one-dimensional 
spin- 1 system showing exclusively ferromagnetic exchange. The spin lattice of CdVOs entails zigzag 
chains with an effective intrachain coupling J ~ —90 K and interchain couplings of J c ~ —18 K and 
J a — —3 K. Quantum fluctuations are partially suppressed by the sizable interchain coupling J c 
that leads to an intermediate regime between one-dimensional and two-dimensional ferromagnetic 
systems. Apart from the peculiar spin model, CdVOs features an unusual mechanism of ferromag- 
netic superexchange. The couplings largely originate from Cd 5s states mediating hoppings between 
half-filled and empty 3d states of V +4 . 

PACS numbers: 75.30.Et, 75.10. Jm, 75.50.Gg, 71.20.Ps 



I. INTRODUCTION 

Ferromagnetism and antiferromagnetism, the two op- 
ponent magnetic interactions, are rarely balanced, be- 
cause specific mechanisms of the magnetic exchange tend 
to favor one of the two options. For example, itiner- 
ant systems are prone to ferromagnetic Stoner instabil- 
ities, whereas superexchange in magnetic insulators is a 
source of mostly antiferromagnetic (AFM) interactions P 
This makes AFM ground states more common among 
insulating transition-metal compounds. Antiferromag- 
netism utterly dominates in low-dimensional magnets, 
where long-range and typically AFM couplings between 
the low-dimensional units (chains or layers) induce the 
overall AFM order. In particular, most of the quasi-one- 
dimensional (ID) spin-| ferromagnetic-chain systems are 
antiferromagnetically ordered because the interchain cou- 
plings are AFMP™ Low-dimensional spin-^ magnets 
with ferromagnetic (FM) ground state are st ill r are and 
restricted to systems based on organic radicals^' or Cu +2 
compounds with non-trivial orbital orderingPHU 

The aforementioned trend is violated by one peculiar 
compound, the low-pressure modification 10 of CdVOs. 
Its ID crystal structure (Fig. [IJp"^ featuring zigzag chains 
of edge-sharing VO5 pyramids, seemingly represents an 
archetypal AFM insulator, where orbital degrees of free- 
dom are eliminated by the square-pyramidal coordina- 
tion of V +4 . The interchain couplings are long-range 
and likely AFM because of the underlying V-O-O-V su- 
perexchange pathways that typically lead to AFM cou- 
plings, while the FM superexchange is usually operative 
at short distances (see, e.g., Refs. IT2lfT4|) . Moreover, the 
next-nearest-neighbor intrachain coupling J 2 between the 
corner-sharing vanadium pyramids should also be AFM, 
as in CaV2C>5 and related compounds (see Sec. |v| and Ta- 
ble |v| ESEZl However, experimental data disprove these 
empirical arguments, and reveal FM order in CdV03 be- 
low Tc = 24 K. 11 Magnetic susceptibility fitted with an 
expression for the classical FM chain model yields an ef- 
fective intrachain coupling J = —100 KpQ Dai et alP^ 
found Ji ~ —288 K and J2 ~ —90 K from generalized 




FIG. 1. (Color online) Crystal structure of CdV03 featuring 
zigzag chains of VO5 pyramids. The chains are separated by 
Cd atoms (spheres). 



gradient approximation (GGA) band structure calcula- 
tions. They tentatively ascribed the unusual FM J2 to 
the shift of the vanadium atom toward the base of the 
pyramid. However, their estimate of J\ severely exceeds 
the experimental coupling of —100 K and, more impor- 
tantly, does not address the nature of the interchain cou- 
plings, which are the driving force of the puzzling FM 
order in CdVOs. 

In the following, we explore the microscopic mechanism 
of ferromagnetism in CdV03 by extensive band structure 
calculations combined with magnetization measurements 
and quantum Monte-Carlo (QMC) simulations. The ap- 
plication of diverse approaches for the evaluation of ex- 
change couplings, along with the direct comparison to the 
experimental data, leads to a reliable microscopic mag- 
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netic model of CdV03. We show that the FM behavior of 
this compound is peculiar, and relate the ferromagnetism 
to Cd 5s orbitals mediating the FM superexchange. Be- 
low, the methodological part (Sec.lnJ is followed by band 



structure results in S ec 
imental data in Sec. |IV 
and summary in Sec. [V 



In] and an analysis of the exper- 
We conclude with a discussion 



II. METHODS 

The band structure was calculated within the 
framework of DFT using the full-potential local- 
orbital scheme (FPL08 . 50-32) PI We applied the local- 
density-approximation (LDA) with the Perdew-Wang 
parametrizatiorP^ for the exchange-correlation potential. 
Exchange couplings were evaluated via two different pro- 
cedures, a model approach and a supercell approach. In 
the former, the LDA band structure was mapped onto 
a multi-orbital Hubbard model that was further treated 
perturbatively in the strongly correlated limit. In the su- 
percell approach, the correlation effects in the V 3d shell 
were treated in a mean-field fashion by the LSDA+/7 
method. Total energies for a number of collinear spin 
configurations were mapped onto a classical Heisenberg 
model to yield individual exchange couplings. Since 
the supercell approach led to somewhat puzzling results, 
we performed an extensive cross-check using: i) GGA 21 
within FPLO; ii) Vienna ab initio simulation package 
(VASP}22l that performs projected augmented wave cal- 
culations and therefore employs a different basis set! 23 ^ 4 ! 
LDA results were obtained for the orthorhombic crystal- 
lographic unit cell, with a fine k mesh of 1040 points in 
the symmetry-irreducible part of the first Brillouin zone. 
DFT+[/ calculations utilized supercells doubled along b 
or c, with k meshes of 150 — 200 points. 

To evaluate the magnetic susceptibility and Curie tem- 
perature of the proposed spin model, we performed QMC 
simulations using the loop algorithm 2 ^ from the ALPS 
package. 26 ID and two-dimensional (2D) finite lattices 
comprised N = 120 and 512 (32 x 16) sites, respec- 
tively, and ensured the absence of finite-size effects for the 
magnetic susceptibility within the temperature range un- 
der investigation. The Curie temperature was estimated 
from simulations for a three-dimensional (3D) model with 
different lattice sizes (see Sec. |IV[). 

Experimental magnetization data were collected with 
a SQUID magnetometer (Quantum Design MPMS) in 
the temperature range 2 — 380 K in applied fields up to 
5 T. The single-phase polycrystalline sample of CdV03 
was prepared by a solid-state reaction of CdO and VO2 
in an evacuated silica tube at 700 °C for 24 hours. A 
25 % excess of CdO was introduced to compensate for 
the losses caused by the volatalization and reaction with 
the tube. The phase purity of the sample was checked by 
x-ray powder diffraction (Hubcr G670 Guinier camera, 
CuK Ql radiation, image-plate detector, 29 = 3 — 100° 
angle range). 
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FIG. 2. (Color online) LDA density of states for CdV0 3 . 
The Fermi level is at zero energy. 



III. BAND STRUCTURE AND EXCHANGE 
COUPLINGS 

A. LDA and model approach 

The LDA band structure of CdV0 3 (Fig. [2]), with oxy- 
gen 2p valence bands below —3 eV and vanadium 3d 
bands at t he Fe rmi level (Ep), is reminiscent of other 
V +4 oxides.E^U The contribution of cadmium is, how- 
ever, large r than typical for an alkaline-eart h (e.g. , Ca +2 , 
Sr +2jEEa or even a d w ( Ag +i^ Zn +2pMm cation. 

The bands below —8 eV originate from the filled Cd Ad or- 
bitals, whereas the states at 3 — 4 eV show predominantly 
Cd 5s character. It is worth noting that the contributions 
of Cd and O at the Fermi level are comparable (6.2% and 
9.6%, respectively, for E < 0.2 eV), yet oxygen 2p states 
dominate over Cd 5s between 1 eV and 3 eV. The sizable 
contribution at Ep distinguishes Cd from other cations 
that also produce conduction bands 3 — 4 eV above the 
Fermi level, but show a negligible contribution at Ep [for 
instance, Pb+ 2 in PbZnVO(P0 4 ) 2 (Ref. HHD or Se+ 4 in 
VOSe03 (Ref. 30)]. The obtained gapless energy spec- 
trum originates from the underestimation of correlation 
effects in LDA. LSDA+f/ reproduces a band gap of about 
2.0 eV (Ud = 4 eV, FPLO) in reasonable agreement with 
the brown color of CdVOs. 

The band complex between —0.5 eV and 4.2 eV com- 
prises 24 bands (Fig. [3]). Since there are four formula 
units per cell, these bands arise from six orbitals per for- 
mula unit: five V 3d and one Cd 5s. The 3d levels of 
vanadium lie below the Cd states, and show a crystal -field 
splitting characteristic of V +4 05 square pyramids! 31 * 32 ! 
Directing the z axis along the short (apical) V-0 bond 
of the VO5 pyramid, we find the lowest-lying d xy crystal- 
field level, represented by four narrow bands at the Fermi 
level (Fig. [3}. 

To analyze exchange couplings, we fit 20 vanadium 
bandd^l with a tight-binding model based on Wannier 
functions adapted to specific orbital symmetries,^ and 
map the resulting transfer integrals t (hoppings)2S! onto 
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FIG. 3. (Color online) Left: LDA band structure for CdV0 3 
(thin light lines) and the fit of the tight-binding model (thick 
dark lines) for the V 3d-related bands. Right: orbital-resolved 
DOS for V 3d states. The Fermi level is at zero energy. The 
poor fit in the vicinity of the T point is related to the strong 
hybridization with the Cd 5s states represented by the four 
high-lying bands (not fitted). The notation of k points is: 
r(0,0,0), X(0.5,0,0), 5(0.5,0.5,0), Y(0, 0.5,0), Z(0, 0,0.5), 
and [7(0.5,0,0.5). 



a multi-orbital Hubbard model with the effective on-site 
Coulomb re pulsio n U e s — 4 eV and Hund's coupling 
J cS = 1 eVP^EH In the t < U cS limit, this model is 
reduced to the Kugel-Khomskii m odel, and the exchange 
couplings are expressed as follows J2Z122I 



U 2 



E 



At 2 



, a Jeff 



([/ off + A Q )(£/ cff + A Q - J cff )' 



(1) 



where t xy and t xy 



are transfers between the 



xy 



states and from the xy (half-filled) to a (empty) states, 
respectively.^ A Q stands for the crystal-field splitting 
between the xy and a states. Note that the states with a 
certain 3d orbital character and the resulting hoppings in 
Eq. ([I ) refer to Wannier functions centered on vanadium 
sites. Each Wannier function entails one of the vanadium 
3d orbitals along with oxygen 2p and cadmium 5s states 
(Fig. [7]). Therefore, metal-ligand transfers are implicitly 
contained in the hopping parameters tj. 

With Eq. Q, one evaluates full exchange couplings 
that are a sum of the AFM superexchange, arising from 
the transfers between the half-filled xy states [Eq. (Ji}, 
first term], and the FM superexchange due to the hop- 
pings to empty d states [Eq. 0, second term]. The 
efficiency of the LDA-based model approach has been 
demonstrated in Refs. [Ml and ESI providing a direct 
comparison to the experiment. 

The leading exchange couplings calculated according 
to Eq. ([I]) are listed in Table [T] The FM components of 
Ji and J2 surpass the AFM superexchange. The inter- 
chain coupling along c lacks any AFM component, yet 
there is a sizable FM contribution of —17 K. The inter- 
chain coupling along a is weaker but also FM. Further 
couplings are below 2 K (in terms of the absolute value) 



TABLE I. Interatomic distances (in A) in the CdVC>3 struc- 
ture and the exchange couplings (in K) calculated with 
Eq. Q: the AFM (J AFM ) and FM (J FM ) contributions and 
the resulting total exchange (J). 



Distance 


jAFM 


jFM 


./ 


Ji 3.05 


18 


-87 


-69 


J 2 3.60 


21 


-30 


-9 


J 3 5.93 





-8 


-8 


J c 5.20 





-17 


-17 


Ja 5.79 


1 


-4 


-3 


with an exception of J3 


~ -8 


K, which is 


the third- 



neighbor intrachain coupling. 

The exclusively FM couplings in CdVC-3 readily lead 
to the FM long-range order. Therefore, the results of 
the model approach are consistent with the experimental 
data, at least on the qualitative level. In Sec. |IV[ we 
will further demonstrate a good quantitative agreement 
with the experiment, while the rest of the present section 
is focused on the application of the supercell approach 



(Sec. Ill B and III C ) and on the microscopic origin of 
ferromagnctism in CdV03 (Sec. 



HID) 



B. DFT+C7 puzzles 

The supercell approach to the evaluation of ex- 
change couplings is based on LSDA+C7 (more gener- 
ally, DFT+C7) calculations that effectively reproduce 
the gapped energy spectrum of a Mott insulator. The 
DFT+J7 method rests upon the mean- field solution of 
the Hubbard model in the strongly correlated limit, but 
the incorporation of this solution into the self-consistent 
procedure leads to several features making the DFT+£/ 
results different from that of the perturbative treat- 
ment on top of LDA (model approach). One differ- 
ence is the application of the on-site Coulomb repul- 
sion and exchange to individual atomic-like orbitals in 
DFT+[/ instead of hybridized LDA bands (correspond- 
ing to molecular-like orbitals or Wannier functions) in 
the model approach. Therefore, the DFT+J7 parameters 
Ud and Jd are g eneral ly different from U e s and J e g in the 
Hubbard modelPffl The second feature is the doublc- 
counting-correction (DCC) scheme that subtracts part of 
the Coulomb energy already contained in LSDA, and en- 
ables the self-consistent procedure. Depending on the fill- 
ing of individual d orbitals, the DCC is applied in either 
around-mean-field (AMF) or fully-localized-limit (FLL) 
fashions. 

The choices of Ud, Jd, and the DCC are made on em- 
pirical grounds, and retain a certain ambiguity. Here, 
we focus on the previously overlooked effect of DCC, 
with Ud and Jd fixed at 3 eV and 1 eV, respectively. 
These values have been justified by FPLO L SDA+L/ calcu- 
lations for several V +4 compounds! 12 * 31 ^ We also varied 
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TABLE II. Exchange couplings (in K) evaluated by the 
DFT+U supercell procedure for different functionals, double- 
counting-correction (DCC) schemes, and band structure 
codes. The DFT+J7 parameters are set to Ud = 3 eV (FPLO), 
U d = 4 eV (VASP), and J d = 1 eV (both codes). 



Jl + J 3 


■h 


Ja 


Jc 


Functional 


DCC 


Code 


-17 


21 


2 


9 


LSDA+(7 


AMF 


FPLO 


-18 


30 


2 


8 


GGA+U 


AMF 


FPLO 


-130 


11 


-4 


-24 


LSDA+ry 


FLL 


FPLO 


-117 


17 


-3 


-17 


GGA+U 


FLL 


FPLO 


-122 


-4 


-3 


-26 


LSDA+C/ 


FLL 


VASP 



Ud in the physically reasonable range of 2.5 — 6 eV, but 
no qualitative differences were found. To avoid calcula- 
tions for large supercells, we estimated J\ + J3 instead of 
evaluating J\ and J3 separately. According to Table [TJ 
|Ja|<C I Ji\, i.e., one may assume J\ + J3 ~ J\. 

The AMF and FLL results for the exch ang e couplings 
in CdV0 3 arc rather different (Table |n|i f23 AMF eval- 
uates the mostly antiferromagnetic scenario, while FLL 
renders J\ strongly FM. We further checked these re- 
sults against the different exchange-correlation potential 
(GGA+U) and the different band structure code. For a 
given DCC, the LSDA+C/ and GGA+U estimates closely 
match. The VASP calculations can be done for FLL 
only, and support the respective estimates from FPLO. 
A marginal difference between the FPLO and VASP re- 
sults (especially for the second-neighbor coupling J2) is 
related to different basis sets and, consequently, differ- 
ent projection schemes employed in the construction of 
the DFT+J7 occupation matrix. This difference also ex- 
plains the 1 eV offset in the value of Ud (see Ref. I59"|) . 
Another remark regards the GGA (without U) results by 
Dai et a^P^who found both J\ and J2 FM, with a clearly 
overestimated absolute value of J\ . Their results cannot 
be directly compared to ours, because uncorrelated GGA 
calculations heavily underestimate correlation effects (for 
example, the reported band gap of about 0.5 e\^is much 
too small to explain the brown color of CdVOs). The ne- 
glect of strong correlation effects typically leads to huge 
errors in the exchange couplings, as demonstrated, e.g., 
in Ref. HQ 

Table [Til shows that the DFT+[7 estimates are robust 
with respect to the exchange-correlation potential and to 
the particular basis set employed in the band structure 
code. Thus, the problem stems from the choice of the 
DCC, which has a strong and unanticipated effect on the 
computed exchange couplings. A simple qualitative anal- 
ysis identified the correct ground state for the FLL set of 
exchange parameters, whereas AMF predicts the wrong 
ground state. Indeed, the FM ground state in CdV03 re- 
quires FM interchain couplings J a and J c and the FM or 
weakly AFM J 2 . The AFM next-nearest-neighbor intra- 
chain coupling J2 frustrates J±, but does not break the 



FM ground state for J 2 /\Jx\ < \ (Ref.HU). Th e FLL re- 
sults J a , J c < and J2/I Ji| < 0.15 fulfill both conditions. 
By contrast, ^2 / 1 ^1 1 — 1 obtained in AMF induces a spi- 
ral order along the chains that are further coupled anti- 
ferromagnetically. Therefore, the AMF-based scenario is 
unrealistic. 

The above analysis puts forward the advantages of 
FLL in evaluating the exchange couplings for CdV03. 
This conclusion is reinforced by the physical meaning 
of different DCC schemes.^ While AMF describes the 
regime of moderate correlations, FLL should be appro- 
priate for strongly correlated systems (^ <C U e s), such 
as CdV03. To check whether this DCC can be used 
universally, we calculated exchange couplings for several 
simple V +4 compounds and further explored the effect of 
the DCC on the stability of different magnetic states in 
DFT+U. 



C. Double-counting correction: AMF vs. FLL 

To compare the performance of the AMF and FLL ver- 
sions of DFT+C7, we select five representative V +4 com- 
pounds showing chain-like magnetic behavior ( Table [ill] ). 
While details of the interchain couplings and the mag- 
netic ground state may be different (and, in some cases, 
not fully understood), the AFM nature of th e intr achain 
couplings is safely established experimentallyP^ESHZlThe 
key difference between our test systems is the mutual ar- 
rangement of the vanadium polyhedra. To analyze the 
connectivity, we only consider the magnetic unit, a VO5 
square pyramid, despite the actual crystallographic local 
environment is often described as an octahedron. The ad- 
vantage of our description is the simple identification of 
the relevant superexchange pathway through the oxygen 
atoms lying in the basal plane of the pyramid. Since all 
systems under consideration reveal the d xy orbital ground 
state, the axial oxygen atom neither connects the pyra- 
mids nor contributes to the superexchange, and only the 
oxygen atoms in the basal plane take part in the superex- 
change couplings. 

The test systems represent several different regimes 
(Fig. [3: i) edge-sharing VO5 pyramids (V0S b204 and 
the low-pressure modification of MgV03 ) j 44 l 45 l h) corner- 
sharing V0 5 pyramids with either twisted (Ba 2 V 3 09, V- 
O-V angle of 96.4° }Por nearly linear (CsVOF 3 , V-O-V 
angle of 164.9° geometries; iii) VO5 pyramids bridged 
by a non-magnetic acetate (CHsCOO - ) groupP^ 

The intrachain exchange coupling J is evaluated as 
the energy difference between the FM and AFM spin 
configurations. Further couplings, both interchain and 
long-range intrachain, are rather weak and irrelevant for 
the present analysis. Quantitative estimates of these 
weak couplings and detailed structural information can 
be found in the preceding studies, where the exchange 
integrals were successfully evaluated using the model 
approach™™ 



In Table III we compare AMF and FLL re- 





Ba 2 V 3 9 CsVOF 3 

FIG. 4. (Color online) Superexchange pathways in the test 
compounds listed in Table |iTl| 

TABLE III. Test quasi- ID V +4 compounds, the V-V dis- 
tances d (in A), the connections between the VO5 pyramids, 
and the exchange couplings J (in K) obtained from the AMF 
or FLL supercell calculations and from the experiment. 



Compound 


d(V-V) 


Connection 


jAMF 


jF LL 


J cxp Ref 


MgVOa 


2.98 


edge-sharing 


128 


-84 


100 




VOSb 2 4 


3.01 


edge-sharing 


248 


154 


250 


ma 


Ba 2 V 3 9 


3.01 


corner-sharing 


79 


84 


94 


EH 


CsVOF 3 


3.91 


corner-sharing 


143 


157 


132 




VO(CH 3 COO) 2 


3.48 


via ac bridges 


544 


528 


430 





suits for the leading intrachain exchange couplings.^ 
The corner-sharing (Ba 2 V30g, CsVOF 3 ) or indirect 
[VO(CH 3 COO) 2 ] connections between the V0 5 pyramids 
lead to a remarkably small difference between AMF and 
FLL. By contrast, there are large discrepancies for the 
couplings between edge-sharing pyramids, especially in 
MgVC>3. These discrepancies are basically independent 
of the Ud value, because J AMF and J FLL show similar 
evolution with a nearly constant offset of J AMF — J FLL ~ 
100 K in VOSb 2 4 and about 200 K in MgV0 3 (Fig. |HJ. 
In VOSb204, the discrepancy might be still tolerable, 
because both types of DCC yield the correct antiferro- 
magnetic solution. The weaker exchange coupling, along 
with a larger offset between AMF and FLL, render the 
FLL result for MgVC»3 qualitatively wrong. Fig. [5] shows 
that any reasonable value of Ud leads to negative J FLL 
violating the experimental AFM coupling. The VASP cal- 
culations for MgVC>3 produce similar J values and con- 
firm the intrinsic nature of the problem. 

Based on our findings for the test compounds, we arrive 
at two conclusions that are - at this stage - rather empiri- 
cal: i) the particular choice of the DCC scheme is relevant 
for short-range couplings only; more specifically, only the 
couplings between the edge-sharing vanadium pyramids 
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FIG. 5. (Color online) Intrachain exchange coupling 
in MgVOa calculated by LSDA+t/ for different values of 
Coulomb repulsion parameter Ud and different DCC schemes: 
AMF (open symbols) and FLL (filled symbols). The marginal 
difference between the FLL results obtained in FPL0 and in 
VASP is likely related to the different basis sets. 

are affected; ii) for edge-sharing pyramids, neither DCC 
can be used universally, because FLL produces a realis- 
tic scenario for CdVC>3, while failing for MgVC>3, and in 
AMF the situation is exactly the opposite. We believe 
that the broad range of coupling geometries considered 
in our study makes this empirical conclusions a helpful 
guidance for future computational work on V +4 oxides 
and other transition-metal compounds. We also empha- 
size the remarkably good performance of the model ap- 
proach for CdVC>3. The model approach is free from the 
double-counting problem, lacks any ambiguity, and rep- 
resents an appealing alternative to the DFT+C7 supercell 
calculations. A further discussion of methodological as- 
pects and tentative remarks on the possible origin of the 
DCC effects are given in Sec. |Vj 

D. Origin of ferromagnetism 

The FM behavior of CdV0 3 contrasts with the AFM 
properties of other V +4 compounds, such as quasi-lD 
MgV0 3 and VOSb 2 4 or quasi-2D CaV 2 5 , CaV 3 7 , 
and CaV40 9 . While J\ is a nearly 90° and possibly FM 
V-O-V superexchange, the weakness of J 2 as well as the 
FM nature of J a and J c are less clear from a microscopic 
point of view. Dai et a/ff^ claimed that the shift of the 
vanadium atom toward the basal plane of the VO5 pyra- 
mid could cause FM J 2 . However, they did not verify 
this conjecture, and refrained from any analysis of the 
FM interchain couplings, which are crucial for the FM 
long-range order. Here, we argue that the origin of fer- 
romagnetism in CdVC>3 is different, and relates to the 
admixture of Cd states to the magnetic orbitals (Wan- 
nier functions). We justify this mechanism by consider- 
ing a model system, a hypothetical CaV03 compound 
with the crystal structure of CdVC>3 but Ca occupying 
the Cd position. In our analysis, we follow the approach 
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FIG. 6. (Color online) LDA density of states for the fictitious 
CdV03-type CaVG-3 (compare to Fig. |5J. The Fermi level is 
at zero energy. 



TABLE IV. Exchange couplings (in K) in the hypothet- 
ical CaVOs compound with the CdVC>3 structure. The 
J AFM , J FM , and J = J AFM + J FM values are obtained from 
Eq. |TJ (model approach), whereas the numbers in the last 
two columns (AMF and FLL) are calculated via the super- 
cell approach with different DCC (LSDA+J7, U d = 3 eV, 
J d = 1 eV). Similar to Table [IT] the LSDA+[7 estimates of 
Ji are in fact Ji + J3 with J3 ~ 5 K according to the model 
analysis. 







Model approach 


LSDA- 


¥U 




jAFM 


J FM 


J 


J 


J 










AMF 


FLL 


Jl 


26 


-42 


-16 


-28 


-163 


Ji 


105 


-23 


82 


133 


83 


Jc 











1 


-3 


Ja 


1 


-1 





2 


-2 



of Refs. QH 123 and ED that utilize the ability of DFT to 
evaluate microscopic parameters of fictitious compounds 
and, therefore, investigate the influence of different struc- 
tural features on the magnetic properties. 

The band structures of CdVC>3 and hypothetical 
CaVC>3 are rather similar, yet CaVC>3 is free from the 
low-lying Cd 5s bands (compare Figs. [2] and |6j. There- 
fore, vanadium bands are largely mixing with O 2p and 
basically lack the cation contribution. 35 This change in 
the band structure has a strong effect on individual 
transfer integrals and the resulting exchange couplings 
(Table [jV| . Compared to Table |IJ we find: i) the in- 
crease in AFM J2; ii) the AFM nature of J a and J c ; 
iii) isotropic interchain couplings (J a ~ J c ). The hypo- 
thetical CaV03 is predicted to be predominantly antifer- 
romagnetic, and conforms to the trends established for 
V +4 compounds. Specifically, the long-range interchain 
couplings are AFM, whereas J2 is about 100 K, as in 
MgV 2 5 and CaV 4 9 (see Sec. [V] and Table [V}. 

The LSDA+f/ results for CaVC>3 reveal the same 





CdVO, 



FIG. 7. (Color online) Vanadium d xy -based Wannier func- 
tions for CdV0 3 (left) and hypothetical CaV0 3 (right). Note 
the Cd 5s contributions (marked with arrows) that are miss- 
ing in the Ca-containing compound. 



strong dependence on the DCC, as previously observed 
in CdV0 3 (compare Tables P and [IV|) . The AMF and 
FLL calculations basically agree on the sizable AFM J2, 
but strikingly differ in the estimate of the short-range 
coupling J\. While we can not obtain any experimental 
information on the hypothetical CaV03 compound, it is 
instructive to compare the LSDA+[7 estimates of J\ to 
the independent result from the model approach. The 
model estimate of Ji ~ — 16 K is now in good agree- 
ment with the AMF prediction Ji ~ —28 K, while the 
FLL estimate of J\ ~ — 163 K is far too large in terms of 
the absolute value. A further discussion of methodologi- 
cal problems related to the DCC of DFT+J7 is given in 
Sec. El 

The role of the Cd 5s states in the magnetism of 
CdV03 is elucidated by Wannier functions. In Fig.[7J we 
compare the V d^-based Wannier functions for CdVC>3 
and CaV03. Apart from the leading V d xy contributions, 
both Wannier functions involve oxygen p orbitals. In 
CdV03, there is an additional Cd 5s contribution, which 
is missing in CaVC>3 . Similar features are found for Wan- 
nier functions based on empty vanadium orbitals. The 
Cd orbitals represent "tails" of the Wannier functions 
and amplify interorbital hoppings that drive the FM su- 
perexchange. 



IV. ANALYSIS OF THE EXPERIMENTAL 
DATA 

The magnetic susceptibility (x) of CdV03 steeply in- 
creases from 380 K to 2 K and clearly indicates the FM 
nature of the system (Fig. [8|. At low temperatures, the 
magnetization reaches the saturated value of about 1 fig 
at the applied field of less than 0.2 T. The lack of hys- 
teresis may be related to a very weak anisotropy of V +4 
(Ref. [52]). Above 230 K, the susceptibility follows the 
Curie- Weiss law 
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FIG. 8. (Color online) Magnetic susceptibility of CdVC>3 
measured in an applied field of 0.5 T (filled circles) and the 
fits with the classical ID model (dashed line) as well as the 
quantum 2D Ji — J c model (solid line). Inset: magnetization 
curve at 2 K. 



where xo — —70 x 10~ 6 emu/mol is the temperature- 
independent contribution of core diamagnetism and van 
Vleck paramagnetism, C = 0.368 emu K mol -1 is the 
Curie constant, and 9 = —46 K is the Weiss tempera- 
ture. The C value corresponds to the effective magnetic 
moment of 1.70 /is that perfectly matches the expected 
value of gnByS(S + 1) = 1.697 ^b with g = 1.96 from 
ESRPEQ The negative Weiss temperature is a signature 
of FM couplings leading to a positive deviation from the 
Curie- Weiss behavior below 230 K. 

Onoda and NishiguchPil fitted the magnetic suscepti- 
bility with the expression for the classical spin chainpS 
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1 the classical AFM chain shows the 
I J ~ 0.35, whereas for 



Our data can be described in a similar way (xo — 
-70 x 10~ 6 emu/mol, g = 1.99, J = -90 K, dashed 
line in Fig. [8]), but the model itself does not apply to 
the spin-| compound CdV0 3 because of inherent quan- 
tum fluctuations in low-dimensional spin-^ systems. For 
example, at S = \ 
susceptibility maximum at T, 
the quantum AFM chain T max /J ~ 0.64. A pronounced 
difference should also be expected for the FM case. In- 
deed, the simulated curve for the quantum FM chain 
(J = — 120 K) fits our experimental data down to 100 K 
only (short-dashed line in Fig. [9]). At lower temperatures, 
the ID quantum model underestimates the susceptibility. 
This should be understood as an effect of quantum fluc- 
tuations that disturb the parallel alignment of spins and 
therefore reduce x- 

Since the experimental data for CdVOa conform to 
the classical model (see Fig. [8]), quantum fluctuations 
in this compound are less pronounced than in the single 
quantum spin chain. This reduction could be related to 
the interchain coupling J c that increases the dimension- 
ality, or the FM intrachain coupling J 2 that increases 
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FIG. 9. (Color online) Fits of the magnetic susceptibility 
with different ID and 2D spin models, see text for details. 



the number of bonds at a lattice site without chang- 
ing the dimensionality. We start with the first option, 
which is also favored by DFT since |J 2 | <C \Ji\- The 
2D Ji — J c model fits the experimental data down to 
30 K with J 1 = -90 K, J c jJ x = 0.2, g = 1.96, and 
Xo = —80 x 10~ 6 emu/mol (solid line in Fig. |9j. For 
comparison, we also considered the isotropic 2D model 
(FM square lattice, J c = J% = —50 K) where quan- 
tum fluctuations are further suppressed by the increased 
dimensionality. This model consequently overestimates 
the susceptibility below 90 K (dash-dotted line in the 
same figure). Therefore, CdVOs exhibits an intermedi- 
ate regime between ID and 2D FM systems. 

The long-range FM order in CdVOa is stabilized by 
the interchain couplings J a and J c . To evaluate the Curie 
temperature Tq , we considered the 3D Jj — J c — J a spin 
model and calculated the Binder ratio of magnetization 
B(T) — (m 4 ) / (m 2 ) 2 for finite lattices of different size. 
The Curie temperature was determined as the crossing 
point of several B{T) curves calculated for L x L/2 x L/2 
finite lattices with L < 64. Similar to Refs.l54landl55"l we 
reduced the dimensions of the lattices along J a and J c to 
account for the anisotropic nature of our system. Using 
J a /J 1 = 0.03 (Tables I and |n]), we arrive at T c fJ\ = 
0.212 (Tq = 19 K) that is slightly below the experimental 
value of 24 KP 

The marginal underestimate of Tq in the J\ — J a — J c 
model may be related to magnetic anisotropy, which lies 
beyond the scope of the present study, or the second- 
neighbor intrachain coupling J2. According to DFT re- 
sults, this coupling is either weakly FM (Table |l| or 
weakly AFM (Table [n). The AFM J 2 introduces frus- 
tration and disables the QMC techniques because of the 
sign problem. By contrast, the effect of FM J2 can be 
readily evaluated. At weak FM J 2 , the overall energy of 
the intrachain exchange (—90 K) is merely redistributed 
between Ji and J 2 . Using J\ = J2 = —40 K, we are 
able to introduce more significant changes by reducing 
quantum fluctuations and improving the susceptibility 
fit even within the purely ID J\ — J 2 model (long-dashed 



line in Fig. |9|. The second- neighbor coupling J2 brings 
the system closer to the classical regime, thereby reduc- 
ing the fitted interchain coupling. Since the ordering 
temperature is mostly sensitive to the value of J c , the 
J 1 — J2 — Ja, — Jc model further underestimates Tc ■ There- 
fore, the sizable FM J2 is unlikely. In conclusion, we ar- 
gue that the unfrustrated J\ — J a — J c spin lattice, provid- 
ing a remarkable fit of the susceptibility and a reasonable 
estimate of Tc, is a valid microscopic model of CdVC>3. 
This model strongly supports our computational results: 
compare the experimental Ji ~ —90 K, J a ~ —3 K, and 
J c ~ -18 K to the calculated J x ~ -69 K, J a ~ -3 K, 
and J c ~ —17 K (Table 0. 



V. DISCUSSION AND SUMMARY 

The peculiar ferromagnetism of CdVC>3 puts forward 
several important issues. First, we have demonstrated 
the importance of the DCC scheme as one of the deli- 
cate parameters of the DFT+J7 method. While it is usu- 
ally difficult to make a well-justified choice of the DCC 
scheme, empirical recipes can be used, as we show be- 
low. We found that FLL produces accurate estimates for 
certain compounds, such as CdV03, but it may fail in 
closely related systems (e.g., MgVOs) where AMF is, on 
the contrary, the method-of-choice. The most pragmatic 
and safe solution to this problem would be the adjust- 
ment of the DCC against the experimental data for each 
specific compound. While this does not diminish the cru- 
cial role of band structure calculations in the microscopic 
modeling of complex materials, the constant reference to 
the experimental data deprives computational methods 
of one important advantage, the ability to predict the 
magnetism of hitherto unexplored systems. To remedy 
this drawback, we performed a test study of several V +4 
compounds with different coupling geometries. 

Based on the comparative study in Sec. |III C| we con- 
clude that the ambiguity related to the choice of DCC is 
confined to short-range couplings only. Moreover, the un- 
certainty is only present for edge-sharing VO5 pyramids 
in MgV0 3 , CdV0 3 , and VOSb 2 4 . The similarly short 
V-V distance of about 3.0 A in P^VsOg with the corner- 
sharing pyramids does not cause any problems: both 
AMF and FLL results are in excellent agreement with 
the experiment (Table |Hl| ) . The principal difference be- 
tween the edge-sharing and corner-sharing geometries is 
the combination of the direct V-V exchange and V-O-V 
superexchange in the former, while only the V-O-V su- 
perexchange is featured by the latter. It is the combi- 
nation of the direct exchange and superexchange or the 
direct exchange itself that are not well reproduced by 
DFT+[/, presumably, due to the oversimplified mean- 
field treatment of correlation effects. The accurate re- 
sults of the LDA-based model approach (Table |l| ensure 
the high accuracy of the hopping parameters used as an 
input. Therefore, the underlying interatomic interactions 
are well reproduced even in LDA, but a better treatment 



TABLE V. Exchange couplings between corner-sharing VO5 
pyramids: V-O-V angles (in deg), V-V distances (in A), and 
exchange couplings J (in K) estimated from the experiment 
(marked with an asterisk) or DFT calculations. 



Compound 


a(V-O-V) 


d(V-V) 


J 


Ref. 


BaoV^Oq 


96.4 


3.01 


94* 


[271 


MgV 2 O s 


117.6 


3.37 


92 






141.1 


3.69 


144 


ED 


CaV 4 9 


129.9 


3.54 


148 


ED 


CaV 2 5 


132.9 


3.49 


608 






135.3 


3.60 


122 


ED 


CdVOs 


136.1 


3.60 


< 20 


This work 


PbV0 3 


147.7 


3.80 


203* 


M 


CsVOF 3 


164.9 


3.91 


132* 


M 



of the on-site correlation effects is required. The devel- 
opment of the respective computational techniques will 
be a rewarding but challenging task that lies far beyond 
the scope of our study. 

Presently, we are able to formulate the following 
recipes for calculating exchange couplings in transition- 
metal compounds. The conventional supercell approach 
(DFT+U calculations for different spin configurations) 
can be safely used for long-range couplings and even for 
short-range couplings not involving the direct overlap of 
the magnetic orbitals. In this case, the choice of the 
DCC plays a minor role compared to other factors, such 
as the adjustment of the Coulomb repulsion parameter 
U4. Once the direct overlap of the magnetic orbitals is 
encountered, DFT+C/ results should be taken with cau- 
tion and tested against experimental data or against in- 
dependent computational estimates. Particularly, we put 
forward the model approach as an appealing alternative 
to the DFT+J7 calculations. The perturbative treatment 
of the multi-orbital Hubbard model is based on the re- 
liable LDA input, free from ambiguity, and uncovers in- 
dividual hopping processes that are responsible for the 
superexchange. We believe that the combination of the 
model and supercell approaches is a viable and reliable 
tool for studying magnetic systems, as demonstrated by 
our present work on CdV03. 

After commenting on the methodological aspects, we 
discuss the physics of CdV03. The unexpected FM be- 
havior of this compound originates from Cd 5s states 
mediating the hoppings between the half-filled and empty 
states of V +4 . The interaction is, therefore, a superex- 
change, but its mechanism is d iffer ent from the conven- 
tional orbital ordering scenario The orbital ordering 
induces magnetic (half-filled) orbitals of different symme- 
try on the neighboring atoms, and favors the hoppings 
between the half-filled and empty orbitals of the same 
symmetry. In CdV03, there is only one type of the mag- 
netic orbital, hence ferromagnetism is driven by hoppings 
between orbitals of different symmetry. 
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The effect of Cd is easily recognized in Table [V] that 
compares exchange couplings between the corner-sharing 
VO5 pyramids. Although one generally expects the in- 
crease in the AFM superexchange at V-O-V angles close 
to 180°, this trend does not hold for V +4 compounds due 
to the relevance of other structural features. For exam- 
ple, basal planes of the pyramids coincide in CaX-^Os and 
CSVOF3, but remain nearly perpendicular in Ba2VaOg 
(Fig. [4| . This explains the broad range of possible V-O- 
V angles and their weak effect on the superexchange. The 
only general trend is the AFM nature of the coupling be- 
tween the corner-sharing VO5 pyramids. The couplings 
typically range between 100 and 200 K, and always ex- 
ceed 90 K. CdVC>3 is a remarkable exception, with the 
AFM coupling reduced below 20 K because of the Cd 5s 
orbitals altering the superexchange. 

The cation-mediated superexchange in vanadium ox- 
ides is not restricted to CdV03. For example, the unusu- 
ally strong and frustrating second-neighbor exchange in 
PbV03 is likely driven by Pb 6p orb itals that marginally 
contribute to the V d xy bandsP^EZl A similar "diagonal" 
superexchange has been proposed for Pbo.55Cdo.4sV205, 
although the intrinsic disorder of Pb and Cd atoms in this 
compound hampered an accurate experimental estimate 
of the respective couplingPHl The aforementioned systems 
are, however, different from CdVC>3, because the Pb and 
Cd cations mediate an AFM superexchange. The non- 
magnetic cations may have diverse effects on the superex- 
change, and the specific scenario is actually determined 
by the hybridization between the cation orbitals and oxy- 
gen orbitals entering vanadium-based Wannier functions. 
In many systems (e.g., vanadium phosphates^), there is 
no effect at all because the cation states are expelled from 
the Fermi level. The cation-mediated superexchange is 
not ubiquitous and requires specific coupling geometries, 
but it can be expected for a variety of cations featuring 
empty s (e.g., Zn+ 2 , Cu +1 , Ag +1 ) or p (e.g., Bi+ 3 , Sn+ 2 ) 
states closely above the Fermi level (see also Ref. I5T))) . 
Respective compounds are likely to host non-trivial spin 



lattices and unusual magnetism. 

Finally, we note that CdV03 is interesting on its own 
as a system showing an intermediate regime between 
the ID and 2D ferromagnets. Experimental studies of 
FM uniform chains helped to find appropriate theoret- 
ical tools for solving the Heisenberg mod el in ID, and 
disclosed peculiar soliton-type excitations! 60 * 61 ! The re- 
spective systems are close to the ID limit, whereas the 
opposite, 2D limit is realized in K2CUF4 and related 
Cu+ 2 halides with layered perovskite-type structures!^! 
CdV0 3 could be a reference point between these two 
qualitatively different model regimes that represent the 
ground states with zero (ID) and non-zero (2D) ordered 
magnetic moment. The crossover between the ID and 2D 
regimes is rather well studied for the AFM casef^ yet a 
comparative study of the FM case could be insightful. 

In summary, we have developed a microscopic mag- 
netic model of CdV03 , and found out the origin of ferro- 
magnetism in this compound. We argue that CdVOs is 
a ferromagnetic spin chain system with an effective (and 
mostly nearest- neighbor) intrachain coupling of —90 K 
and a sizable interchain coupling J c = — 18 K along one 
of the dimensions. Our model is based on extensive band 
structure calculations and verified by a direct comparison 
to the experimental magnetic susceptibility and Curie 
temperature. The unusual ferromagnetic couplings arise 
from Cd 5s orbitals that contribute to the vanadium- 
based Wannier functions, mediate hoppings between the 
half-filled and empty d states of vanadium, and lead to 
the ferromagnetic superexchange. This mechanism puts 
forward diverse effects of non-magnetic cations on su- 
perexchange in transition-metal compounds. 
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